UNCLASSIFIED 


AD  401  342 


(lepAaduced, 
lufr  the. 

DEFENSE  DOCUMENTATION  CENTER 


FOR 


SCIENTIFIC  AND  TECHNICAL  INFORMATION 


CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  sped* 
flcatlons  or  other  data  are  used  for  any  purpose 
other  than  In  connection  vlth  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  whatsoever]  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  right b 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


■Jb  ’ 


oi 


i 


401  342 


OF  THE  TOTAL  THERMAL  RADIANT  ENERGY 
^"*"^|lIrED  BY  AN  UNDERWATER  EXPLODING  WIRE 

s 

J***' 

cz 


CD  j.  8.  Hege 


J.  8- 

c 5 


o 


GO 


63-3^ 

DASA  13M» 
U8NRDL-TR-612 
10  January  1963 


1 


U.S.  NAVAL  RADIOLOGICAL 
DEFENSE  LABORATORY 

SAN  FRANCISCO  77.  CALIFORNIA 


1 2ND.  P7483 

I 


1 


RADIOLOGICAL  EFFECTS  BRANCH 
E.  A.  Schuert,  Head 


CHEMICAL  TECHNOLOGY  DIVISION 
L*  H-  Gevantaan,  Head 


AMINI8TRATIVE  INFORMATION 


The  work  reported  was  done  as  part  of  a ;  pro¬ 
ject  sponsored  by  Defense  Atomic  Support  Agency , 
under  NWER  Program  A7a,  Subtask  10.001.  The 
project  is  described  (as  Program  A-l,  Problem  2) 
in  USNRDL  Technical  Progam  Summary  for 
Fiscal  Years  1963 ,  1964.  and  1965.  1  November 
1962. 


Eegene  P.  Cooper 
Scientific  Director 


^2$ 

E.  8.  Rofh,  CAPT  USN 
Commanding  Officer  and  Director 


ABSTRACT 


The  total  thermal  radiation  from  an  underwater  spark  has  been 
determined  by  the  measurement  of  the  absolute  radiation  at  two  wave¬ 
lengths  as  a  function  of  time.  Planck's  law  for  thermal  radiation  was 
applied. 

The  spark  was  generated  by  electrically  exploding  a  5-mil.  copper 
wire  stretched  between  two  submerged  electrodes.  The  spark  was  assumed 
to  be  cylindrical  with  its  length  determined  by  the  length  of  the  wire. 

The  light  measurements  were  made  with  two  calibrated  photoelectric 
tubes.  Wavelengths  of  407  and  6l0  millimicrons  were  selected  with 
Farrand  interference  filters. 

Values  for  the  rates  of  total  thermal  radiation  and  values  for 
the  spark  temperatures  and  radii  were  tabulated  and  graphed  as  func¬ 
tions  of  time  for  three  explosions  under  three  different  sets  of  con¬ 
ditions.  Total  thermal  radiation  was  obtained  by  graphically  integrat¬ 
ing  the  rate  of  thermal  radiation  over  the  time  of  the  process.  For 
two  of  the  three  explosions  approximately  30  percent  of  the  spark  energy 
was  emitted  as  thermal  radiation.  It 'was  not  possible  to  determine 
this  fraction  for  the  third  explosion. 

Analysis  of  the  propagation  of  error  shows  that  the  error  in  the 
results  for  total  thermal  radiation  is  less  them  7*5  times  any  error 
in  the  measurements  of  the  radiation  incident  on  the  photoelectric  tubes. 


SUMMARY 


Problem 


The  energy  of  an  underwater  spark  le  dissipated  in  the  shock  wave, 
in  thermal  radiation  (light  energy)  and  thermal  conduction,  and  in  the 
steam  bubble.  The  energy  from  an  underwater  nuclear  explosion  is  dis¬ 
sipated  via  the  same  pathways.  The  energy  partition  among  these  path¬ 
ways  is  uncertain.  If  measurements  of  the  total  thermal  radiation 
emitted  by  an  underwater  spark  could  be  made,  the  values  obtained  may 
then  be  useful  in  experiments  with  miniature  models  of  underwater 
nuclear  explosions.  They  could  also  be  compared  with  the  estimates  of 
the  radiant  energy  emission  from  actual  underwater  nuclear  explosions. 

Findings 

Values  for  the  spark  temperature,  its  radius  and  its  rate  of  radi¬ 
ant  energy  emission  were  found  as  functions  of  time.  From  these  the 
total  thermal  radiation  was  determined.  It  was  found  that  about  3°  per¬ 
cent  of  the  spark  energy  was  emitted  as  thermal  radiation. 
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INTRODUCTION 


To  investigate  the  possibilities  of  simulating  certain  aspects  of 
an  underwater  nuclear  explosion  with  a  miniature  explosion,  NRDL  is 
studying  some  of  the  processes  occurring  when  a  wire  is  electrically 
exploded  underwater.  A  major  problem  in  the  study  of  the  exploding 
wire  is  to  determine  the  explosion  energy  partition  among  the  result¬ 
ant  shock  wave,  steam  bubble  and  electromagnetic  radiation.  This  report 
is  devoted  specifically  to  the  determination  of  the  total  thermal  radi¬ 
ant  energy  emitted  during  the  brief,  brilliant  flash  of  the  explosion. 

The  experimental  approach  to  this  problem  is  the  conventional  one 
where  radiation  from  the  explosion  is  measured  at  two  wavelengths  as  a 
function  of  time  with  calibrated  photoelectric  tubes.  A  value  for  the 
total  thermal  radiant  energy  can  then  be  computed  using  Planck's  law 
for  radiation  from  a  black  body  source.  Two  questions  arise  at  this 
point:  (l)  Is  the  system  radiating  with  a  Planckian  distribution? 

(2)  Is  the  method  of  analysis  so  sensitive  to  errors  in  the  measurements 
that  irreducible  experimental  error  causes  a  prohibitively  large  error 
in  the  final  solution? 

Substantial  evidence  Is  given  by  Martin  (Bef.  l)  in  favor  of  the 
system's  radiating  as  a  black  body.  He  demonstrated  theoretically  that 
the  underwater  spark  can  be  considered  as  progressing  continuously 
through  successive  states  of  Thermal  equilibrium  so  that  photons  inci¬ 
dent  on  this  spark  surface  are  distributed  in  Planckian  manner.  Martin 
further  offered  experimental  evidence  that  the  spark  surface  radiates 
like  a  black  body  in  the  visible  region  of  the  spectrum.  Since  light 
measurements  in  the  present  experiment  are  made  in  the  visible  region  of 
the  spectrum  it  seems  justifiable,  at  least  for  the  purposes  of  this 
experiment,  to  treat  the  spark  as  a  black  body. 

The  second  question  is  investigated  in  Section  IV.  It  is  shown 
that  the  error  in  the  value  obtained  for  the  total  radiant  energy  emit¬ 
ted  is  less  than  7*5  times  any  error  incurred  in  the  measurement  of  the 
experimental  parameters.  It  is  felt  that  this  is  surprisingly  small 
and  that  the  approach  outlined  above  is  feasible. 
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SWITCH 


Pig.  1  Schematic  Diagram  of  the  Wire -exploding  Facility 


The  underwater  spark  is  produced  in  an  indoor,  water-filled  tank, 
with  a  copper  wire  stretched  "between  two  centrally  located  electrodes 
(Fig.  l).  Electrical  energy  is  switched  from  a  charged  capacitor  bank 
into  the  circuit  containing  the  wire.  The  charged  capacitor  and  attend- 
ent  circuitry  which  includes  the  wire  is  an  LCR  system.  When  the  switch 
is  closed,  the  resultant  current  is  a  damped  sinusoidal  wave  (Fig.  2) . 
The  wire  is  vaporized  and  a  conducting  plasma  of  copper  atoms  and  water 
molecule  fragments  is  formed.  The  shape  of  this  plasma  has  been  shown 
by  photographic  methods  to  be  cylindrical  with  the  axis  of  the  cylinder 
defined  by  the  wire. 2  The  plasma  is  called  the  spark  channel  and  the 
radiation  is  emitted  from  the  spark  channel  surface.  The  resistance  of 
the  spark  channel  comprises  approximately  40  $  of  the  total  resistance 
of  the  circuit  so  that  about  40  #  of  the  energy  originally  stored  in 
the  capacitor  is  dissipated  through  the  spark  channel. 3 


METHOD  OF  ANALYSIS  AMD  WORKING  EQUATIONS 


The  detailed  method  of  analysis  used  to  determine  the  total  thermal 
radiant  energy  emitted  at  the  spark  channel  surface  is  given  below. 

First  the  determination  of  the  channel  temperature  is  discussed  both 
from  theoretical  and  experimental  points  of  view.  The  Stefan  -Bo It zmann 
law  is  then  introduced  to  give  the  rate  of  energy  radiation  per  unit 
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Fig-  2  Oscillograph  of  a  Typical  Channel  Current 
Gap  =  1  inch;  stored  energy  =  1800  j 
Oscilloscope  vertical  sensitivity  =  1  v/cm 
Oscilloscope  horizontal  sensitivity  =  20  nsec /cm 


area  of  channel  surface.  The  method  of  determining  the  total  channel 
surface  area  i6  treated,  and  the  equations  for  the  rate  of  thermal 
radiation  from  the  entire  spark  channel  surface  is  derived-  Finally, 
the  working  equations  are  presented  by  substituting  known  or  measured 
values  for  the  constants.  The  total  thermal  energy  radiated  during 
the  whole  process  is  obtained  upon  graphical  integration  of  the  rate 
of  energy  emission. 

Spark  Channel  Temperature 

The  determination  of  the  temperature  of  a  black  body  from  values 
for  the  radiation  emitted  at  two  different  wavelengths  is  the  basis 
of  the  experimental  method  used.  Planck's  law  for  black  body  radiation 
requires  that  the  radiation  emitted  at  wavelength  \  per  unit  time, 
per  unit  surface  area,  per  unit  wavelength,  be 

he 

t  o  .2,-5  /  XkT  .x-1 

=  2  jt  he  \  \e  -  1) 
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where  h  =  Planck's  constant 
c  =  speed  of  light 
X  =  wavelength  of  the  radiation 
k  =  Boltzmann's  constant 

T  =  black  body  temperature  in  degrees  Kelvin 


The  ratio  of  two  J^'s  defines  a  unique  temperature  for  a  black  body. 


If  the  temperature  of  a  black  body  is  varying  in  time,  then  meas¬ 
urement  of  the  J^'s  for  two  wavelengths  as  functions  of  time  defines 
the  time  function  of  the  temperature.  It  is  by  these  means  that  the 
time  function  of  the  spark  channel  temperature  is  determined. 


The  ratio  /j^  is  denoted  by  the  symbol  B.  Solving  the  equation 


B  = 


j  _5  /  — \-l 

\  V  le\kT  -1 1 


J>  X'5 


he 
Xr&T 


e  -  1 

for  T  gives  the  spark  channel  temperature. 


The  actual  measurement  of  the  radiation  at  a  given  wavelength  is 
made  by  a  photoelectric  tube  which  has  a  Farrand  interference  filter 
covering  the  cathode  to  absorb  all  light  except  the  narrow  band  at  the 
desired  wavelength.  The  photoelectric  tube  is  calibrated  so  that  the 
absolute  value  of  the  radiant  flux  incident  on  the  cathode  can  be  read 
directly  from  the  voltage  generated  by  the  tube.  The  voltages  are  re¬ 
corded  on  an  oscillograph  over  the  time  of  the  emission  process. 


The  relationship  between  J\,  the  radiation  emitted  at  the  spark 
channel  surface,  and  E^,  the  radiant  flux  incident  on  the  photoelectric 
tube  cathode,  can  be  written  as  follows: 

2  y  x 

x  D2 

where  R  is  the  spark  channel  radius,  L  is  the  channel  length,  and  D  is 
the  distance  from  the  spark  channel  to  the  photoelectric  tube.  The  tube 
must  be  aligned  so  that  it  sights  along  the  perpendicular  bisector  of 
the  spark  channel  axis  (Fig.  3)  and  D  must  be  much  greater  then  R. 

The  factor  2  RL  in  the  numerator  is  the  cross-sectional  area  of  the 
channel  which  is  seen  by  the  tube  cathode.  The  symbol  7.  is  the 
attenuation  factor  of  the  media  between  the  spark  channel  and  the  photo¬ 
electric  tube  for  radiation  ai;  wavelength  X. 
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Fig.  3  Arrangement  for  Photoelectric  Tubes  and  Spark  Channel 


The  value  of  7^  is  alvays  a  number  less  than  one  except  when  a 
perfect  vacuum  exists  between  the  channel  and  the  tube.  It  is  assumed 
that  the  attenuation  in  air  is  insignificant  for  the  distance  involved 
so  that  7^  need  be  determined  only  for  the  water  and  glass  in  the  explo¬ 
sion  tank.  A  value  for  7\  is  obtained  by  measuring  the  flux  of  wave¬ 
length  k  incident  on  the  photoelectric  tube  from  a  light  source  placed 
at  the  position  of  the  wire  in  the  tank,  first  with  the  tank  empty  and 
the  glass  window  open,  and  then  with  the  tank  full  and  the  window  closed. 
The  ratio  of  the  latter  measurement  to  the  former  is  7\. 

B,  the  number  which  defines  the  spark  channel  temperature ,  can  now 
be  found  in  terms  of  the  experimental  parameters,  and  7\,  as  follows: 


This  expression  is  true  only  when  both  photoelectric  tubes  (one  for 
measuring  the  radiation  at  each  wavelength)  are  properly  aligned  at  the 
same  distance  from  the  spark  channel.  Then, 
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Total  Thermal  Radiant  Energy 

Once  T  is  computed  from  B,  the  rate  at  which  energy  is  radiated 
per  unit  area  from  the  channel  surface  for  all  wavelengths  is  given  by 
the  Stefan -Boltzmann  law-  Where 

4 

W  =0-  T 


and 

0~  =  the  Stefan -Boltzmann  constant 

W  multiplied  by  the  channel  surface  area,  A^,  gives  P,  the  rate  of 
energy  radiation  for  the  entire  channel  surface. 

Ac  can  be  found  from  the  equation, 

2  “  \  ^ 


in  the  following  manner: 


A  =  2  it  RL 
c 


Jv  7  \ 

*  \  =  (2  jt  RL)  - g- 

rt  D 


2  2 


Substituting 


for  gives 


2  *  hcV5  (e1^1  -  I)’1 
n  D2  \5  hc/XkT 

Ac =  I : '  r  (e 


1) 
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Then  P,  the  rate  at  which  energy  is  radiated  from  the  entire  spark 
channel  surface  is: 


P 


<r-  « 


(ehc/^.1)T\ 

2  he 


The  total  radiant  energy  emitted  by  the  channel  is  the  integral  of  P  dt 
taken  over  the  entire  explosion  process.  This  integral  is  evaluated 
graphically  by  computing  enough  values  for  P  to  describe  a  curve  of  P 
versus  time  and  then  determining  the  area  under  the  curve  with  the 
planimeter. 

The  Working  Equations 

E|y  using  predetermined  values  for  the  calibration  constants  and 
the  load  resistances  of  the  photoelectric  tubes,  and  by  substituting 
numbers  for  the  physical  constants,  the  equations  can  be  written  in  the 
working  form  actually  used  in  reducing  the  experimental  data  to  useful 
information. 


The  current  i^  developed  in  the  photoelectric  tube  due  to  radiant 
flux  incident  on  the  cathode  is  related  to  K  by  the  equation  i^  = 
where 7^  is  the  calibration  constant.  The  voltage  V^,  across 
tne  load  resistance,  R^,  is  R^ix  and  is  recorded  on  the  oscillograph. 
The  ratio,  B,  then,  can  be  written  in  terms  of  the  voltages  recorded  as 


B 


V 


v'o: 


Using  the  values  for  the  R's  andT^'s  which  are  given  in  the  section 
"Instrumentation." 


B  =  0.424 


where  k.  is  407  mu  and  kg  16  6l0  mp  in  all  cases  in  this  report.  To 
facilitate  the  calculation  of  the  temperature,  T,  from  a  value  for  B, 
a  standard  curve  of  B  versus  T  is  graphed,  using  the  equation 
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c  hc/X  kT 

*|(e  2  -1) 

B  '  >.[  (.‘“Ajte  .  T) 

This  graph  is  shown  in  Pig.  k. 

The  equation  for  the  rate  at  which  energy  is  radiated  from  the 
spark  channel  surface, 

P  -  -r°y  ^  -  1)  T4 

2  he  ^ 

can  be  evaluated  using  the  data  for  either  wavelength,  407  mu  or  6l0  mil. 
In  theory  both  sets  of  data  should  give  the  same  value  for  P.  In  prac¬ 
tice,  however,  the  values  are  not  exactly  the  same  so  that  both  are 
evaluated  and  their  average  value,  P  =  (p^q^  +  ?6io )/®,  *6  usei  to  graph 
the  relationship  P  versus  time .  * 

The  working  equations  for  P^Qy  and  Pg10  sure 

P407  -  °'m  *  ^  \07  (e3^  *  ^-l)  T4  ergs/sec 

P610  -  1-57  *  M"12  j2—  Vg10  (e2,36  x  10VT.l)  i4  ergs/sec 

6io 


where  D,  the  distance  from  the  spark  channel  to  the  photoelectric  tube, 
is  in  centimeters  and  the  voltage  on  the  oscillograph,  is  millivolts. 

The  spark  channel  radius  can  be  computed  from  the  equation 

k 

P  =  A  cr  T 
c 

as  follows: 

Ac  =  2  it  RL  =  — tt 
<r  T 


- E 

2  it  L<rT 
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centimeters 


R  = 


2.8  x  10 

JT 


L  T 


where  P  is  in  ergs  per  second  and  L  is  in  centimeters. 

The  quantities  described  above  are  tabulated  and  graphed  in 
Section  V  for  3  explosions  under  3  different  sets  of  conditions. 


INSTRUMENTATION 


The  photoelectric  tubes  used  for  making  the  radiation  measurements 
at  407  and  6l0  mu  are  the  RCA  929  and.  925  tubes  respectively.  The  407 
and  6l0  mil  Farrand  interference  filters  have  peak  transmissions  of  31 
percent  and  23  percent  and  half  bandwidths  of  17  mu  and  13  mu  respectively. 
In  addition  to  the  interference  filter,  a  neutral  filter  is  used  with 
the  929  photoelectric  tube  in  order  to  keep  the  tube  current  below  levels 
at  which  saturation  might  occur.  The  neutral  filter  has  a  transmission 
of  approximately  5  percent  at  407  mp. 

Figure  5  is  a  schematic  diagram  of  the  photoelectric  tube  circuits. 
Each  tube  is  powered  by  a  3°0  volt  dry  cell  battery  which  has  a  0.101  pf 
capacitor  across  its  terminals  to  insure  against  a  momentaiy  loss  of 
voltage  in  the  case  of  large  currents  in  the  photoelectric  tube  circuits. 
Load  resistances  of  532  ohms,  2.21  kilohms,  21.4  kilohms,  215  kilohms 
and  2  megaohms  are  available  in  the  circuit  for  the  tube  measuring  radi¬ 
ation  at  407  mp.  For  the  other,  load  resistances  of  535  ohms,  2.23  kil¬ 
ohms,  20.5  kilohms,  210  kilohms, and  2  megohms  can  be  selected.  The 
voltages  generated  across  the  load  resistances  Eire  measured  on  a  single 
Tektronix  dual  beam  type  551  oscilloscope  with  two  type  K  fast -rise -time 
preamplifiers.  During  calibration  and  during  the  determination  of  the 
7\'s,  the  more  sensitive  type  D  pre -amplifiers  are  used,  because  the 
radiation  intensities  are  much  lower. 

Each  photoelectric  tube  with  its  circuit  components  is  shielded 
electrostatically  in  an  aluminum  box  in  which  the  window  is  covered 
by  the  filter.  The  boxes  are  carefully  sealed  against  all  light  except 
that  entering  through  the  filter. 
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TEMPERATURE  CKxlO4) 


Fig.  h  Ratio  B  of  Flux  Intensities  at  Wavelengths  707  and  6l0  mp. 

from  a  Black  Body 


Fig.  5  Schematic  Diagram  of  Photoelectric  Tube  Circuitry 
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Each  filter -photoelectric  tube  unit  is  calibrated  to  give  radiant 
flux  in  direct  proportion  to  the  current  flowing  through  the  load  resis 
tance.  Once  calibration  is  accomplished  the  only  component  which  can 
be  altered  in  the  unit  is  the  load  resistance. 


Two  independent  methods  of  calibration  are  used  and  the  results 
averaged.  One  method  employs  a  standard  .lamp  which  emits  light  with  a 
flux  density  at  25  cm  from  the  envelope  of  1.6  pw  (cm^.mp) “1  at  407  mp, 
and  13*1  (iw  (cm^'mp)-!  at  6l0  mp  when  powered  by  15.4  amps  A.C.  The 
values  for  flux  densities  are  rms.  The  second  method  employs  a  G.E. 
type  18A-T10-2-6V  lamp  which  has  a  rectangular  tungsten  filament  with 
a  surface  area  on  one  side  of  0.144  cm2.  The  apparent  temperature  of 
the  filament  is  determined  with  an  optical  pyrometer  and  the  true  temp¬ 
erature  is  found  using  values  for  the  spectral  emissivity  (e^)  of  tung¬ 
sten  from  Forsythe  and  Worthing.^  The  flux  density  at  407  and  6l0  mp, 

D  cm  from  the  filament  on  an  axis  normal  to  the  plane  of  the  filament 
at  its  center,  is  calculated  from  the  equation 


«  D2 


where  Af  is  the  filament  area  and  where  J}.  is  computed  using  Planck's 
law.  The  value  for  the  true  temperature,  and  the  values  for  are 
taken  from  Forsythe  and  Worthing.  Since  the  flux  density  and  the  photo 
electric  tube  current  are  related  by  the  equation  ty  =7^] 2^,^,  the 
calibration  constant  for  the  filter -photoelectric  tube  unit  measuring 
the  radiation  at  wavelength  X,  can  be  evaluated  from  the  equation 


'I 


X  = 


«  D2  i. 


\  JX  Af 


The  two  methods  agreed  within  7  percent  and  gave  average  values  of 

=  0.088l  pa*cm2«mp*  (p)~\ 

^610  =  °*°392  pa-cm2-mp-(p)_1. 


In  selecting  the  proper  load  resistances,  there  are  two  factors  to 
consider.  One  is  to  have  the  load  resistance  great  enough  so  that  the 
stray  pick-up  detected  by  the  photoelectric  tube-oscilloscope  circuitry 
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that  is  due  to  the  large  currents  in  the  exploding  wire  facility,  is 
not  significant  compared  to  the  voltage  across  the  load  resistance. 

The  other  is  to  have  the  RC  time  constant  (load  resistance  multiplied 
by  the  input  capacitance)  small  enough  to  make  it  certain  that  the 
oscilloscope  is  responding  accurately  to  the  light  signals  incident 
upon  the  photoelectric  tube.  The  Tektronix  type  K  pre -amplifiers  with 
their  low  input  capacitance  of  20  (4tf  are  used  with  only  22  inches  of 
amphenol  type  RG-62A  connecting  cables.  The  tubes  are  placed  close 
enough  to  the  exploding  wire  that  the  photoelectron  currents  are  as 
great  as  possible  without  reaching  saturation  levels.  This  requires 
that  load  resistances  of  about  20  kilohms  be  used.  The  actual  values 
of  the  load  resistances  were  determined  with  a  Leeds  and  Northrup  ohm- 
meter. 

An  attempt  was  made  to  determine  how  large  a  current  the  photo¬ 
electric  tubes  could  carry  for  times  of  100  microseconds  before  their 
responses  became  non-linear.  Measurements  of  the  light  from  identical 
explosions  at  progressively  shorter  distances  from  the  spark  channel 
were  made.  The  currents  recorded  were  up  to  50  microamperes.  No 
attempt  was  made  to  measure  currents  higher  than  50  microamperes  so 
that  the  actual  saturation  currents  may  have  been  considerably  higher 
for  those  short  times.  Currents  measured  during  actual  data  runs  were 
below  these  values. 


THE  PROPAGATION  OP  ERROR 


The  precision  with  which  experimental  determinations  are  known 
depends  upon  two  prime  factors.  First,  the  technical  precision  with 
which  the  experimental  parameters  are  measured,  and  second,  the  manner 
in  which  any  error  in  the  measurement  of  the  parameters  is  propagated 
via  the  mathematics  used  to  develop  the  final  expression. 

The  standard  deviation  of  a  number  of  determinations  is  one  way  of 
expressing  precision,  and  it  embodies  the  two  factors  mentioned  above. 

Since  the  number  of  determinations  in  the  present  experiment  were 
limited  no  standard  deviations  could  be  found;  only  qualitative  state¬ 
ments  can  be  made  about  precision.  These  qualitative  statements  are 
based  upon  the  good  reproducibility  of  the  experimentally  measured  para¬ 
meters  in  the  formative  stages  of  the  experiment,  and  upon  the  determi¬ 
nation  of  an  upper  limit  to  the  growth  of  the  error  in  the  E^'s,  as 
these  parameters  were  mathematically  converted  into  the  expression  for 
total  energy  radiated. 
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The  determination  of  the  upper  limit  to  the  growth  in  error  is 
made  according  to  Scarborough's  method  of  propagation  of  error  and  is 
given  below. 

Scarborough^  shows  that  the  relative  error  in  a  function  having 
several  independent  parameters 

Q  =  Q  (Ux,  Ug,  ...  ug 
is  expressed  by  the  relationship 


E 


A&  _  1 
r  =  Q  Q 


AU, 


+ 


+ 


The  function  of  interest  here  is  the  integral,  /Pdt,  the  total 
thermal  energy  radiated.  The  relative  error  in  thi it  function  can  be 
evaluated  in  terms  of  the  radiant  fluxes,  E407  ^  E6l0>  attenua¬ 
tion  factors,  7407  811(1  ?6l0>  811(1  the  distance,  D,  from  the  spark  channel 
to  the  photoelectric  tubes.  Thus, 

Q  =  J&lt  =  Q  74q7’  76l0' 


and 


_  1  /  & 
Q  \  ^ 


407 


AE 


407 


AD 


It  is  clear  that  the  relative  error  A/pdt/ /Pdt  is  always  equal  to  or 
less  than  the  maximum  value  of  AP/P  *uurirtg  the  time  of  integration  so 
that  the  maximum  value  of  AP/P  provides  an  upper  limit  to  A/pdt//Pdt- 
It  is  this  function  (AP/P)  max.  that  will  be  evaluated.  J  J 


Let 


P  p  ^E4o7'  E6io>  7407>  76io’ 


errors  in  E4Q7  and  Eg10  are  assumed  to  be  the  only  ones  of  significance. 
Then 
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ll 

Since  P  =  a~ AqT  where  Ap  is  spark  channel  surface  area  and  T  is  its 
temperature,  then  by  partial  differentiation 

AP  =  +  4  AT 

P  =  A  T 
c 

Expressing  AAC/AC  and  AT/T  in  terms  of  AE^q^/E^q^  and  ^E6io/E6lO  8iveB 
AP/P  as  a  function  of  the  significant  relative  errors  in  the  experimen¬ 
tal  parameters. 

Evaluation  of  AT/T 

Figure  4  gives  the  relationship  between  the  temperature ,  T,  and 
the  experimentally  determined  number  B.  Any  small  variation  in  T  is 
related  to  a  variation  in  B  by  the  equation 


AB 


—  A 
dT 


where  dB/dT  is  the  slope  of  the  curve  B  versus  T  in  the  region  of  the 
temperature  variation.  Multiplying  by  l/T(dB/dT)-^  gives 


But 


AT  1  (dBl 
T  "  T  l  dT/ 


AB. 


y6l0  E407 
7407  E6io 


and  the  relative  error  in  B  is 


AB  ^407  _  ^610 
B  e407  e6io 

Assuming  that  positive  and  negative  AE^  's  sure  equally  likely  and  that 

^407  ^610 

E407  E6l0 

then 

AB  a*  E^ 

B  S  — 
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and 


AB 


vhere  either  E^q~, 
Substituting 


or  malr  be  used  for  E^. 
this  for  AB  in  the  equation 


for  AT/T  gives. 


AT  2B  dB  "1  ^ 
T  S  T  dT 


Both  B  and  dB/dT  can  be  found  using  the  graph  of  B  versus  T  (Fig.  4). 
Evaluation  of  AAC/AC 

It  has  been  shovn  that 


A  = 
c 


„  D2  ^  j. 

A 


-i) 


2  he2  7. 


The  relative  error  in  A^.  due  to  errors  in  and  T  is  also  derived 
according  to  the  method  of  Scarborough,  and  is 


ehc/\KT  \  AT 
^7^® mlJ  T 

where  Al/T  is  as  above. 

The  expressions  for  AT/T  and  AAC/AC  are  now  substituted  into  the 
original  equation 


AA  AE, 


he 

XKT 


4  AT 
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to  yield 


AP 

P 


he  hc/kkT 
XkT  e 

ehc/\kT  _  ~ 


AEX 

&T 


A  discussion  of  the  evaluation  of  AP/P  can  be  found 


in  the  next  section. 


RESULTS 

Total  Thermal  Radiation 


The  radiation  from  three  exploding  wires  was  measured  by  the  photo¬ 
electric  tubes  and  recorded  on  oscillographs,  and  values  for  the  spark 
channel  temperatures,  radiant  powers,  and  radii  were  computed  as  func¬ 
tions  of  time.  The  times  of  radiation  emission  were  about  100  micro¬ 
seconds.  The  total  thermal  energies  radiated  were  determined  graphically. 
The  conditions  are  tabulated  in  Table  1.  The  computed  quantities  for 


TABIE  1 

Explosion  Conditions 


Explosion 

V* 

,  o. 
(kv) 

E** 

(joule) 

L 

(cm) 

D 

(cm) 

7407 

76io 

1 

25 

1800 

3-81 

730 

0.594 

0.539 

2 

25 

1300 

2.54 

730 

0.563 

0.500 

3 

15 

648 

2.54 

457 

0.526 

0.561 

*  V0  =  Initial  voltage  on  capacitor. 

**£0  =  Energy  initially  stored  in  capacitor. 


all  three  explosions  are  listed  in  Table  2.  The  oscillographs  taken  for 
the  three  explosions  are  shown  in  Fig.  6.  Figures  7,  8,  and  9  are  graphs 
of  the  computed  quantities  over  the  times  of  the  explosions. 
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Computed  Quantities  From  Three  Explosions 


a.  Explosion  1 

Upper  trace  sensitivity  =0-5  V/cm 
Lower  trace  sensitivity  =  0.1  v/cm 
Sweep  speed  =  10  psec/cm 


Bb.  Explosion  2 

Upper  trace  sensitivity  =0-5  V/cm 
Lower  trace  sensitivity  =  0.05  V/cm 
Sweep  speed  =  10  psec/cm 


c.  Explosion  3 

Upper  trace  sensitivity  =  0.2  V/cm 
Lower  trace  sensitivity  =  0.05  v/cm 
Sweep  speed  =  10  psec/cm 


Fig.  6  Oscillographs  of  the  Radiation  Intensities 
Upper  trace  measures  radiation  at  Uo?  mp 
Lower  trace  measures  radiation  at  610  mp 
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Fig.  7  Rate  of  Energy  Radiation  from  the  Spark  Channel  Surface 

for  Three  Explosions 


Fig.  8  Spark  Channel  Temperature  for  Three  Explosions 
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Fig.  9  Spark  Channel  Radius  for  Three  Explosions 


Two  of  the  explosions  were  made  with  the  5.76-nf  capacitor  initially 
charged  to  25  kv,  one  with  a  gap  (wire  length)  of  I-1/2  inches  and  the 
other  with  a  gap  of  1  inch.  The  energy  initially  stored  in  the  capacitor 
in  these  two  cases  was  1800  joules.  For  the  third  explosion,  the  capa¬ 
citor  was  charged  to  15  kv  and  the  gap  was  1  inch.  The  initially  stored 
energy  was  &8  joules. 

In  all  cases  the  wire  was  5  mil  copper,  the  depth  of  the  gap  below 
the  surface  of  the  water  was  6  inches,  the  pressure  over  the  water  was 
ambient,  and  the  water  temperature  was  approximately  75°F* 

Evaluation  of  Error 


The  results  show  that  the  majority  of  the  radiation  occurs  when  T 
is  between  4000  and  15,000°K.  In  this  region  AP/P  is  a  maximum  when  T 
is  in  the  neighborhood  of  15,000°K.  Evaluation  AP/P  when  k  =  6l0  mu 
gives 


AP 
—  < 

P  a 


0-75 


AE610 

E6l0 


The  relative  error  in  the  total  energy  radiated  is  then  less  than  7*5 
times  any  error  in  the  measurements  of  the  radiant  fluxes.  Further  it 
is  estimated  that  these  errors  in  measurement  are  less  than  10  percent. 
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DISCUSSION 


The  total  radiant  energies  emitted  were  15*0,  11.6  and  5*9^  percent, 
respectively  for  explosions  1,  2,  and  3>  of  the  energies  initially  stored 
in  the  capacitors.  These  percentages  are  derived  from  the  values  listed 
in  Tables  1  and  2.  Buntzen,  with  the  apparatus  used  in  the  present 
experiment,  determined  that  42  percent  of  the  initially  stored  energy  is 
put  into  the  spark  channel  when  the  conditions  of  explosion  2  exist. 

Thus,  the  energy  thermally  radiated  in  explosion  2  was  11.6/0.42  or  27-6 
percent  of  the  total  energy  dissipated  by  the  spark  channel,  the  remain¬ 
ing  72.4  percent  going  into  the  shock  wave,  the  bubble  and  into  thermal 
conduction . 

Explosion  1  differed  from  explosion  2  only  in  that  the  spark  channel 
was  1-1/2  times  longer.  Assuming  that  the  resistance  of  the  spark  channel 
varies  directly  with  its  length,  it  can  be  shown  that  if  42  percent  of 
the  stored  energy  went  into  the  spark  channel  in  explosion  2,  then  51 
percent  must  have  gone  into  the  spark  channel  in  explosion  1.  Then  the 
thermal  radiation  was  15/0.51  or  29  percent  of  the  total  spark  channel 
energy.  Fbr  explosions  1  and  2,  roughly  30  #  of  the  energy  put  into  the 
spark  channel  was  later  emitted  as  radiant  energy. 

The  conditions  of  explosion  3  were  identical  to  those  of  explosion  2, 
except  that  in  the  former  the  energy  initially  stored  in  the  capacitors 
was  only  one  third  of  what  it  was  in  the  latter.  The  total  spark  chan¬ 
nel  energy  was  unknown  and  presently  no  explanation  can  be  offered  for 
the  comparatively  low  yield  of  radiant  energy. 

Some  apparent  discrepancies  in  the  experiment  are  indicated  by  the 
values  obtained  for  the  spark  channel  temperatures  and  radii.  These 
show  the  temperatures  increasing  and  the  radii  decreasing  at  times  late 
in  the  emission  process.  It  is  known  that  the  energy  input  to  the  spark 
channel  has  fallen  to  relatively  low  values  at  these  times  (see  Fig.  2) 
so  it  seems  that  the  channel  temperature  should  be  decreasing.  Further, 
there  is  no  reason  for  the  channel  radius  to  decrease  at  a  time  when  the 
pressures  within  are  still  very  high  as  indicated  by  Martin's  work.1 

To  best  evaluate  these  discrepancies  it  is  important  ta. note  how 
error  in  r*  and  error  in  A_  contribute  to  this  error  in&AJIr-  or  P. 

The  final  expression  for  AP/P  in  the  section  "The  Propagation  of  Error" 
contains  the  factor 


21 


k  .(*£_  „W«r)(  „)>=/*»  .  -1 


It  can  be  seen  that  the  number  4  is  from  the  above  expression  for  the 
error  in  T*  while 

j  h£_  ehc/\kTj  j  ehc/\kT  _  ^  -1 


is  from  the  expression  for  the  error  in  Aq.  These  two  quantities  are 
opposite  in  sign,  and,  at  15,000°K,  close  to  being  equal  in  magnitude. 
Therefore  the  relative  error  AP/P  is  less  than  the  relative  error  in 
either  P  or  A,,,  and  the  confidence  with  which  the  values  for  the  total 
radiated  energy  can  be  accepted  is  greater  than  for  either  the  values 
for  Tr  or  for  the  values  for  the  radii  which  are  directly  proportional 
to  Ac. 

Further  investigation  should  be  directed  toward  the  explanation  of 
the  results  for  the  temperatures  and  radii.  Also  the  emissivities  of 
the  spark  showed  surface  should  be  determined  over  as  wide  a  wavelength 
spectrum  as  possible  so  that  it  can  be  known  more  precisely  how  closely 
the  spark  resembles  a  black  body. 
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Commandant,  Command  and  General  Staff  College 

Commandant,  Arny  War  College 

AIR  FORCE 


Deputy  Chief  of  Staff,  Operations  (Operations  Analysis) 
Director  of  Research  and  Development,  DCS/D 
Assistant  Chief  of  Staff,  Intelligence  (AFCIN-3A2A) 
Deputy  Chief  of  Staff,  Operations  (AFOOP) 

Assistant  Chief  of  Staff,  Intelligence,  USAFE 
CG,  Aeronautical  Systems  Division  (ASAPRDONS) 

CG,  Air  Force  Systems  Command  (RDRWA) 

Directorate  of  Civil  Engineering  (AFOCE) 

Director,  USAF  Project  SAND 


27 


1 

1 

1 

1 

1 

2 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


6 

1 

6 

1 

1 

3 

1 

2 
1 
2 
1 
1 
1 
1 
20 
1 
1 
1 
1 


1 

2 

1 

1 

1 

3 

1 

1 


Director,  USAF  Project  RAND  (Dr.  Rapp) 

Commandant,  School  of  Aerospace  Medicine,  Brooks  AFB 
CG,  Strategic  Air  Command  (OAWS) 

Office  of  the  Surgeon  (SUP3.1),  Strategic  Air  Command 
Office  of  the  Surgeon  General 
CG,  Special  Weapons  Center,  Kirtland  AFB 
Director,  Air  University  Library,  Maxwell  AFB 
Commander,  3415th  F&SG 

Commander,  Electronic  Systems  Division  (CRZT) 

Air  Force  Cambridge  Research  Laboratory  (CRQST-2) 

CiC,  Air  Defense  Command  (ADLDC-A) 

Commander,  Tactical  Air  Command 
Commander,  Alaskan  Air  Command  (AAOTN) 

CiC,  Pacific  Air  Forces  (PFCIE-MB) 

Air  Force  Technical  Applications  Center  (Scheid) 
Commander,  Western  Development  Division  (WDSIT) 
Commander,  1009th  Special  Weapons  Squadron 

OTHER  POD  ACTIVITIES 


Chief,  Defense  Atomic  Support  Agency  (Library) 

Chief,  Defense  Atomic  Support  Agency  (Vickery) 

Commander,  Field  Command,  DASA,  Sandia  Base 
Commander,  Field  Command,  DASA,  Sandia  Base  (FCDV) 

Commander,  Field  Command,  DASA,  Sandia  Base  (FCTG) 

Commander,  Field  Command,  DASA,  Sandia  Base  (FCWT) 

Commander,  Field  Command,  DASA,  Sandia  Base  (FCTG5,  Library) 

Office  of  Civil  Defense,  Washington 

Office  of  Civil  Defense,  Washington  (Miller) 

Civil  Defense  Unit,  Army  Library 
Director  of  Defense  Research  and  Engineering 
Director,  Armed  Forces  Institute  of  Pathology 
Commandant,  Armed  Forces  Staff  College 
Commandant,  Industrial  College  of  the  Armed  Forces 
Armed  Services  Technical  Information  Agency 
Director,  Armed  Forces  Radiobiology  Research  Institute 
Director,  Advance  Research  Projects  Agency 
Director,  Weapons  Systems  Evaluation  Group 
Commander,  JTF-8 

AEC  ACTIVITIES  AND  OTHERS 

Research  Analysis  Corporation 

Albuquerque  Operations  Office 

Argonne  National  Laboratory 

AEC  Scientific  Representative,  France 

AEC  Scientific  Representative,  Japan 

Atomic  Energy  Commission,  Washington  (Library) 

Atomic  Energy  Commission,  Washington  (Corsbie,  DBM) 

Atomic  Energy  Commission,  Washington  (Starbird,  DMA) 
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Atomic  Energy  Commission,  Washington  (Hooper) 

Atomic  Energy  Commission,  Washington  (Joseph) 

Atomic  I&iergy  Commission,  Washington  (Hollister) 

Atomic  Energy  of  Canada,  Limited 
Atomics  International 
Battelle  Memorial  Institute 
Beers,  Roland  F.,  Inc. 

Brookhaven  National  Laboratory 

Bureau  of  Mines,  College  Park 

Bureau  of  Mines,  Laramie 

Bureau  of  Mines,  Washington 

California  Institute  of  Technology,  Pasadena 

Carnegie  Institution 

Division  of  Peaceful  Nuclear  Explosives,  Washington 

Division  of  Raw  Materials,  Washington 

duPont  Company,  Aiken 

duPont  Company,  Wilmington 

Fundamental  Methods  Association 

General  Atomic  Division 

General  Electric  Company,  Cincinnati 

General  Electric  Company,  Richland 

General  Electric  Company,  San  Jose 

Grand  Junction  Office 

Harvey  Mudd  College 

Hawaii  Marine  Laboratory 

Holmes  and  Narver,  Inc. 

Home  stake  Mining  Company 

Johns  Hopkins  University  (Wolman) 

Knolls  Atomic  Power  Laboratory 

Lamont  Geophysical  Laboratory 

Los  Alamos  Scientific  Laboratory  (Library) 

Lovelace  Clinic 

Michigan  State  University  (Triffett) 

Moran,  Proctor,  Mueser  and  Rutledge 
Mound  Laboratory 

NASA,  Scientific  and  Technical  Information  Facility 
Nevada  Operations  Office 

Office  of  Assistant  General  Counsel  for  Patents 

Phillips  Petroleum  Company 

Power  Reactor  Development  Company 

Public  Health  Service,  Washington 

Reynolds  Electrical  and  Engineering  Company,  Inc. 

Sandia  Corporation,  Albuquerque  (Library) 

Sandia  Corporation,  Albuquerque  (Merritt) 

Sandia  Corporation,  Albuquerque  (Reed) 

Sandia  Corporation,  Livermore 

San  Francisco  Operations  Office 

Scripps  Institution  of  Oceanography  (Van  Dorn) 

Scripps  Institution  of  Oceanography  (Folsom) 

Space  Technology  Laboratories,  Inc. 

State  of  California,  Department  of  Fish  and  Game  (Bissell) 
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Tennessee  Valley  Authority 

Tracer lab,  Inc.,  Richmond 

Union  Carbide  Nuclear  Company  (ORGDP) 

Union  Carbide  Nuclear  Company  (ORNL) 

U.S.  Coast  and  Geodetic  Survey,  San  Francisco 

U.S.  Coast  and  Geodetic  Survey,  Washington 

U.S.  Geological  Survey,  Albuquerque 

U.S.  Geological  Survey,  Denver 

U.S.  Geological  Survey,  Menlo  Park 

U.S.  Geological  Survey,  Naval  Weapons  Plant 

U.S.  Geological  Survey  (Nolan) 

U.S.  Geological  Survey,  Washington 

U.S.  Geological  Survey,  WR  Division 

U.S.  Mission  to  the  International  Atomic  Qiergy  Agency 

U.S.  Weather  Bureau,  Las  Vegas 

U.S.  Weather  Bureau,  Washington 

University  of  California  (Libby) 

University  of  California,  Berkeley  (Library) 

University  of  California,  Livermore  (Library) 
University  of  California,  Livermore  (Vay  Shelton) 
University  of  California,  Los  Angeles 
University  of  California,  Los  Angeles  (Larson) 
University  of  Hawaii  (Cox) 

University  of  Rochester  (Hempelmann) 

University  of  Tennessee  (UTA) 

University  of  Washington  (Donaldson) 

University  of  Washington  (Flemming) 

Woods  Hole  Oceanographic  Institute 
Yale  University  (Sears) 

Technical  Information  Extension,  Oak  Ridge 

USNRDL 


USNRDL,  Technical  Information  Division 
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